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Abstract 

We show that proper time, when defined in the quantum theory of 2d 
gravity, becomes identical to the stochastic time associated with the stochas- 
tic quantization of space. This observation was first made by Kawai and col- 
laborators in the context of 2d Euclidean quantum gravity, but the relation 
is even simpler and more transparent in he context of 2d gravity formulated 
in the framework of CDT (causal dynamical triangulations) . 
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1 introduction 



Since time plays such a prominant role in ordinary fiat space quantum field 
theory defined by a Hamiltonian, it is of interest to study the role of time 
even in toy models of quantum gravity where the role of time is much more 
enigmatic. The model we will describe in this article is the so-called causal 
dynamical triangulation (CDT) model of quantum gravity. It starts by pro- 
viding an utraviolet regularization in the form of a lattice theory, the lattice 
link length being the (diffeomorphism invariant) UV cut-off. In addition the 
lattice respects causality. It is formulated in the spirit of asymptotic safty, 
where it is assumed that quantum gravity is described entirely by "conven- 
tional" quantum field theory, in this case by approaching a non-trivial fixed 
point [1], [2] . It is formulated in space-times with Lorentzian signature, 
but the regularized space-times which are used in the path integral defining 
the theory allow a rotation to Euchdean space-time. The action used is the 
Regge action for the piecewise linear geometry. Each (piecewise linear) ge- 
ometry used in the path integral has after rotation to Euclidean signature an 
Euclidean Regge action, related to the original Lorentzian action in the same 
way as when one in flat space-time rotates Lorentzian time to Euclidean time 
(see [3, 4] for details of the Regge action in the CDT approach). By rotat- 
ing each Lorentzian CDT lattice (or piecewise linear geometry) to Euclidean 
signature the non-perturbative path integral is performed by summing over 
a set of Euclidean lattices originating from the Lorentzian lattices with a 
causal stucture, and this set is different from the full set of piecewise hnear 
Euclidean lattices. Like in ordinary lattice field theories we approach the 
continuum by fine-tuning the bare coupling constants. The rotation to Eu- 
clidean space-time makes it possible to use Monte Carlo simulations when 
studying the theory, and in four-dimensional space-time, which for obvious 
reasons has our main interest, there exists a region of coupling constant 
space where the infrared behavior of the universe seen by the computer is 
that of (Euclidean) de Sitter space-time [5, 6] (for a pedagodical review, see 
[7]). Non-trivial UV properties have been observed [8], properties which have 
been reproduced by other "field theoretical" approaches to quantum gravity 
[9, 10]. 

Numerical simulations are very useful when trying to understand if a 
non-perturbatively defined quantum field theory has a chance to make sense. 
However, numerical simulations have their limitations in the sense that they 
will never provide a proof of the existence of a theory and it might be difficult 
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in detail to follow the way the continuum limit is approached since it requires 
larger and larger lattices. It is thus of interest and importance to be able to 
study this in detail, even if only in a toy model. Two-dimensional quantum 
gravity is such a toy model which has a surprisingly rich structure. Many of 
the intriguing questions in quantum gravity and in lattice quantum gravity 
are still present in the two-dimensional theory. We will discuss the solution 
to two-dimensional CDT in the rest of this article and we will see that in 
time, which is introduced as proper time, has an interpretation as stochastic 
time in a process there the evolution of space can be viewed as a stochastic. 

2 The CDT formalism 

We start from Lorentzian simplicial spacetimes with d — 2 and insist that 
only causally well-behaved geometries appear in the (regTilarized) Lorentzian 
path integral. A crucial property of our explicit construction is that each 
of the configurations allows for a rotation to Euclidean signature, as men- 
tioned above. We rotate to a Euchdean regime in order to perform the sum 
over geometries (and rotate back again afterward if needed). We stress here 
(again) that although the sum is performed over geometries with Euclidean 
signature, it is different from what one would obtain in a theory of quan- 
tum gravity based ah initio on Euclidean spacetimes. The reason is that not 
all Euchdean geometries with a given topology are included in the "causal" 
sum since in general they have no correspondence to a causal Lorentzian 
geometry. 

We refer to [3] for a detailed description of how to construct the class of 
piecewise linear geometries used in the Lorentzian path integral in higher di- 
mensions. The most important assumption is the existence of a global proper- 
time fohation. This is illustrated in Fig. 1 in the case of two dimensions. We 
have a sum over two-geometries, "stretching" between two "one-geometries" 
separated a proper time t and constructed from two-dimensional building 
blocks. In Fig. 2 we have shown how to fill the two-dimensional space-time 
between the space (with topology S^) at time tn and time tn+i = t^, -\- a 
where a denotes the lattice spacing. While we in the lattice model often 
use units where everj^hing is measured in lattice length (i.e. the lattice links 
have length one), we are of course interested in taking the limit a ^ to 
recover continuum physics. 

In the path integral we will be summing over all possible ways to connect 
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Figure 1: Piecewise linear space-time histories 1+1 dimensional quantum 
gravity 




Figure 2: The propagation of a spatial slice from time t to time t + 1. The 
end of the strip should be joined to form a band with topology x [0, 1]. 

a given Id "triangulation" at time tn and a given Id triangulation at tn+i to 
a slab of 2d space-time as shown in Fig. 2, and in addition we will sum over 
all Id "triangulations" of at times t„. Thus we are demanding that the 
time-slicing is such that the topology of space does not change when space 
"evolves" from time tn to time tn+i- 

The Einstein-Hilbert action S^^ in two dimensions is trivial since there is 
no curvature term as long as the topology of space-time is unchanged (which 
we assume presently). Thus the (Euclidean) action simply consist of the 
cosmological term: 

Sf = \j A\ ^ Sf^'' = (1) 
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where N2 denotes the total number of triangles in the two-dimensional trian- 
gulation. We denote the discretized action the Regge action since it is a trivial 
example of the natural action for piecewise linear geometries introduced by 
Regge [11]. The dimensionless lattice cosmological coupling constant A will 
be related to the continuum cosmological coupling constant A by an additive 
r enormalization : 

A = Ao + ^Aa^ (2) 

the factor 1/2 being conventional. The path integral or partition function 
for the CDT version of quantum gravity is now 

Gf(L„L„T) = J]J-e-^r-(r), (3) 

where the summation is over all causal triangulations T of the kind described 
above with a total of T time steps, an "entrance loop" of length li — Lia and 
an "exit loop" of length I2 = ^20. The factor 1/Cr is a symmetry factor, 
given by the order of the automorphism group of the triangulation T. 

One can, somewhat surprisingly, evaluate the sum over triangulations in 
(3) analytically [12]. It is a counting problem and thus we introduce the 
corresponding generating function. In our model the generating function has 
a direct physical interpretation. We define 

Gf\x,,X2;t) = J2 e-^^^^e-^^^^Gf (Li,L2;T). (4) 

Li,L2 

Thus G^^\Xi, X2; T) is the generating function of the numbers G*^-*(Li, L2; T) 
if we write Zi = e~^^, Z2 = e~^^. But we can also view X as a (bare) 
dimensionless boundary cosmological constant, such that a boundary cosmo- 
logical term X ■ L has been added to the action. In this way G^^\Xi, X2] T) 
represents the sum over triangulations where the lengths of the boundaries 
are allowed to fluctuate, the fluctuations controlled by the values X-i of the 
boundary cosmological constants. In general we expect, just based on stan- 
dard dimensional analysis, the boundary cosmological constants Xi to be sub- 
jected to additive renormalizations when the continuum limit is approached. 
Like (2) we expect 

Xi = Xc + Xi a, (5) 
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where x then denotes the continuum boundary cosmological constant, and 
one, after renormahzation, has the continuum boundary cosmological action 

X • I. 

We refer to [12] for the explicit combinatorial arguments which allow us 
to find (5^^ (Xi, X2; T). Let us just state the following results: one can derive 
an exact iterative equation (using notation Z — e~^, W — , Q — e~^) 

6f(Z,W;T)^^Gf{^,W;T-i) (6) 

This equation can be iterated and the solution found [12]. However, it is easy 
to see that Qc = 1/2 and that Zc = I and we can now take the continuum 
limit in (6) using t — T ■ a and find 



^&;:>ix,y;t) + ^[{x' - \)Gfix,y;t)\ ^0, (7) 

This is a standard first order partial differential equation which should be 
solved with the boundary condition 

(a;,2/;t = 0) = ^ (8) 
X + y 

corresponding to 

G^^\h,l2;t^0)^S{h-l2). (9) 

The solution is thus 

where x{t;x) is the solution to the characteristic equation 

^^-(x'-X), x{t = 0)^x. (11) 
We thus have an explicit solution for Gf\x,y]t) since we obtain 

m = vx (^+")-^"^^^(^-") (12) 

(V^ + a;)+e-2^^*(VA-a;)■ ^ ^ 
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If we interpret the propagator -^(Zi, l2;t) as the matrix element between 
two boundary states of a Hamiltonian evolution in "time" t, 

Gf\h,k;t)=<h\e-^"%> (13) 

we can, after an inverse Laplace transformation, read off the functional form 
of the Hamiltonian operator Hq from (7), 

Hoi^) = ^{^'- A)) , m) = -l^, + A/. (14) 

This end our short review of basic 2d-CDT. We have here emphasized that 
all continuum results can be obtained by explicit solving the lattice model 
and taking the continuum limit simply by letting the lattice spacing a — 0. 
The same will be true for the generalized CDT model described below, but 
to make the presentation more streamlined we will drop the explicit route 
via a lattice and work directly in the continuum. 



3 Generalized CDT 

It is natural the ask what happens if the strict requirement of "classical" 
causality on each geometry appearing in the path integral is relaxed. While 
causality is a resonable requirement as an outcome of a sensible physical the- 
ory, there is no compleling reason to impose it on each individual geometry 
in the path integral, since these are not physical observables. We used it, 
inspired by [13], as a guiding principle for obtaining a path integral which 
is different from the standard Euclidean path integral, which was seemingly 
a necessity in higher than two space-time dimensions since a "purely" Eu- 
clidean higher dimensional path integral did not lead to interesting continuum 
theories. 

In Fig. 3 we show what happens if we allow causality to be violated locally 
by allowing space to split in two at a certain time t, but wc never allow the 
"baby" universe which splits off to come back to the "parent" universe. The 
baby universe thus continues its life and is assumed to vanish, shrink to noth- 
ing, at some later time. We now integrate over all such configurations in the 
path integral. From the point of view of Euclidean space-time we are simply 
integrating over all space-times with the topology of a cylinder. However, 
returning to the original Lorentzian picture it is clear that at the point where 
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Figure 3: In all four graphs, the geodesic distance from the final to the 
initial loop is given by t. Differentiating with respect to t leads to eq. (15). 
Shaded parts of graphs represent the full, (^s-dependent propagator and disc 
amplitude, and non-shaded parts the CDT propagator. 



space splits in two the light-cone is degenerate and one is violating causal- 
ity in the strict local sense that each space-time point should have a future 
and a past light-cone. Similarly, when the baby universe "ends" its time 
evolution the light-cone structure is degenerate. These points thus have a 
diffeomorphism invariant meaning in space-times with Lorentzian structure, 
and it makes sense to associated a coupling constant gg with the process of 
space branching in two disconnected pieces. 
The equation corresponding to Fig. 3 is [14] 

-Gx,gA^,y;t) = -_[((a;2 - A) +2(7, Wx,gM)Gx,gA^,y;t)'\. (15) 

Wx^g^{x) is denoted the disk amplitude with a fixed boundary cosmological 
constant x. It is related to the disk amplitude with a fixed boundary length 

by 

/■oo 

Wx,,Sx)= / dle-^'W,,,M- (16) 
Jo 

It describes the "propagation" of the a spatial universe until it vanishes in 
the vacuum. If we did not allow any spatial branching we would simply have 

Wi'\x)= dtGf\x,l = 0;t) = y=, (17) 

Jo x+vA 

where G^^\x,l;t) denotes the Laplace transform of Gf\l',l;t) with respect 
to From the composition rules for Gx,g^{h,h]t) it follows that it has 
(mass) dimension 1. Thus Gx^g^{x,l2]t) is dimensionless and it follows that 
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Figure 4: Graphical illustration of eq. (58). Shaded parts represent the full disc 
amplitude, unshaded parts the CDT disc amplitude and the CDT propagator. 



the (mass) dimension of the coupling constant Qs must be 3. In a discretized 
theory it will appear as the dimensionless combination ggd^, a being the 
lattice spacing, and one can show that the creation of more than one baby 
universe at a given time t is suppressed by powers of a (see [14] for details). 
Thus we only need to consider the process shown in Fig. 3. For a fixed 
cosmological constant A and boundary cosmological constants x, y expres- 
sions like Gyg^{x, y; t) and W\^g^{x) will have a power series expansion in the 
dimensionless variable 

- (18) 



A3/2 

and the radius of convergence is of order one. Thus the coupling constant 
Qs indeed acts to tame the creation of baby universes and if Qs exceeds this 
critical value eq. (15) breaks down and is replaced by another equation cor- 
responding to Liouville quantum gravity with central change c = (see [14] 
for a detailed discussion). 

Differentiating the integral equation corresponding to Fig. 3 with respect 
to the time t one obtains (15). The disc amplitude Wx^g^{x) is at this stage 
unknown. However, one has the graphical representation for the disc ampli- 
tude shown in Fig. 4. It translates into the equation [14] 

Wx,gM = Wi'\x)+ (19) 

oo oo 

g,j dt j dhdh {li + l2)Gf{xM + l2\t)Wx,gXh)Wx^gXk) 



The superscript (0) indicates the CDT amplitudes without baby universe 
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branching, calculated above. We assume 

W^,,^=,{x) = wf{x), (20) 

and similarly for The integrations in (19) can be performed and if we 

define Wx,gXx) by 

one can show that Wx^g^x) is given by the expression 



= Kx - u)^J{x + uy-2K, (22) 

where 

X = x\f\, - w + K = 0. (23) 

In order to have a physically acceptable Wx^g^x), one has to choose the 
solution to the third-order equation which is closest to 1 and the above 
statements about the expansion of Wx,gXx) in a power series in n follow. 
Eq. (15) can now be written as 

-Gx,gA^,y;t) = -—^Wx,gMGx,gAx,y;t)\. (24) 
In analogy with (7) and (10), this is solved by 

G,,.(.,,;t) = %M-^, (25) 

Wx,gSx) x{t,x)+y 

where x{t,x) is the solution of the characteristic equation for (24), the gen- 
eralization of eq. (11): 

^ = -Wx,gM, x{0,x)^x, (26) 

such that 

dy 



(27) 

m Wx,gM 

This integral can be expressed in terms of elementary functions and one can 
thus find an exphcit expression for Gx,gs{x,y;t) in the same way as eq. (12) 
led to an explicit solution for the Gf\x,y;t) appearing in (10) . 
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4 The matrix model representation 



The formulas (22) and (21) are standard formulas for the resolvent of a 
Hermitean matrix model, calculated to leading order in A^, the size of the 
matrix. In fact the following matrix model 

Z(A, gs) ^jd<t> e-^^^(^), F(0) = ^ - \<t>') (28) 
has a resolvent 

^ (^) ) = ^^'^»(^)(^) + 0{l/N% (29) 

where Wx^g^{x) is given by (21), and where the expectation value of a matrix 
expression C(0) is defined as 

{Om = -^^^ I dct> e-^^^(^) 0{4>). (30) 

What is surpising here, compared to "old" matrix model approaches to 2d 
Euclidean quantum gravity, is that the large N limit reproduces directly the 
continuum theory. No scaling limit has to be taken. The situation is more 
like in the Kontsevich matrix model, which directly describes continuum 
2d gravity aspects. In fact the qubic potential is "almost" like the qubic 
potential in the Kontsevich matrix model, but the wold-sheet interpretation 
is different. 

Can the above correspondance be made systematic in an large N expan- 
sion and can the matrix model representation help us to a non-perturbative 
definition of generalized 2d CDT gravity ? The answer is yes [15]. 

First we have to formulate the CDT model from first principles such that 
we allow for baby universes to join the "parent" universe again, i.e. we have 
to allow for topology changes of the 2d universe, and next we have to check 
if this generalization is correctly captured by the matrix model (28) [16]. 

5 CDT string field theory 

In quantum field theory particles can be created and annihilated if the process 
does not violate any conservation law of the theory. In string field theories 
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one operates in the same way with operators which can create and annihi- 
late strings. From the 2d quantum gravity point of view we thus have a 
third-quantization of gravity: one-dimensional universes can be created and 
destroyed. In [17] such a formalism was developed for non-critical strings 
(or 2d Euclidean quantum gravity). In [16] the formalism was applied to 2d 
CDT gravity leading to a string field theory or third quantization for CDT, 
which allows us in principle to calculate any amplitude involving creation 
and annihilation of universes. 

Let us briefly review this formalism. The starting point is the assumption 
of a vacuum from which universes can be created. We denote this state |0) 
and define creation and annihilation operators: 

[vl>(/),v[/t(/')]=M(/-/'), vI/(/)|0) = (0|vI/t(/) = 0. (31) 

The factor / multiplying the delta-function is introduced for convenience, see 
[16] for a discussion. 

Associated with the spatial universe we have a Hilbert space on the pos- 
itive half-line, and a corresponding scalar product (making Ho{l) defined in 
eq. (14) hermitian) : 

ii^M^ I jm)Mi)- (32) 

The introduction of the operators \I'(/) and ^I'^(/) in (31) can be thought of 
as analogous to the standard second quantization in many-body theory. The 
single particle Hamiltonian Ho defined by (14) becomes in our case the "sin- 
gle universe" Hamiltonian. It has eigenfunctions ■0n(O with corresponding 
eigenvalues e„ = 2ny/X, n — 1,2,...: 

V'n(0 = I e-^'pn-iil), Ho{l)Ml) = e„Vn(0, (33) 

where p„-i(/) is a polynomial of order n — 1. Note that the disk amplitude 
W^\l), which is obtained from (17), formally corresponds to n = in (33): 

Wi'\l) = e-^', Ho{l)wi'\l) = 0. (34) 

This last equation can be viewed as a kind of Wheeler-deWitt equation if 
we view the disk function as the Hartle-Hawking wave function. However, 
W^\l) does not belong to the spectrum of Ho{l) since it is not normalizable 
when one uses the measure (32). 
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Figure 5: Graphical illustration of the various terms in eq. (38). 

We now introduce creation and annihilation operators a}^ and a„ cor- 
responding to these states, acting on the Fock- vacuum |0) and satisfying 
[a„, al^] = 6n,m- We define 

vl>(/) = J^anMl), = Y.^irAl). (35) 

n n 

and from the orthonormality of the eigenfunctions with respect to the mea- 
sure dl/l we recover (31). The "second-quantized" Hamiltonian is 

/■°° r// 

Ho= T ^H0^o(/)^(/), (36) 







.(0) 



and the propag ator G^;;\h,l2;t) is now obtained as 

(/i,/2;t) = (0|vI/(/2)e-*^°vl/t(/^)|o). (37) 

While this is trivial, the advantage of the formalism is that it automati- 
cally takes care of symmetry factors (like in the many-body applications in 
statistical field theory) both when many spatial universes are at play and 
when they are joining and splitting. We can follow [17] and define the fol- 
lowing Hamiltonian, describing the interaction between spatial universes: 

H = Ho -gsjdhj dl2^\h)^\h)^{h + k) (38) 
-ag^j dh J dk^^ih + l2)^ik)^ih)- J jpmil), 

where the different terms of the Hamiltonian are illustrated in Fig. 5. Here 
Qs is the coupling constant we have already encountered in Sec. 3 of mass 
dimension 3. The factor a is just inserted to be able the identify the action 
of the two ^fg-terms in (38) when expanding in powers of Qs- We will think of 
a = 1 unless explicitly stated differently. When a = 1 H is hermitian except 
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for the presence of the tadpole term. It tells us that universes can vanish, but 
not be created from nothing. The meaning of the two interaction terms is as 
follows: the first term replaces a universe of length /1+/2 with two universes of 
length li and ^2- This is one of the processes shown in Fig. 5. The second term 
represents the opposite process where two spatial universes merge into one, 
i.e. the time-reversed picture. The coupling constant gs clearly appears as a 
kind of string coupling constant: one factor Qs for splitting spatial universes, 
one factor Qs for merging spatial universes and thus a factor g'^ when the 
space-time topology changes, but there is also factors for branching alone. 
This is only compatible with an Euclidean SFT-picture if we associate a 
puncture (and thus a topology change) with the vanishing of a baby universe. 
As discussed above this is indeed not unnatural from a Lorentzian point of 
view. From this point of view the appearance of a tadpole term is more 
natural in the CDT framework than in the original Euclidean framework in 
[17]. The tadpole term is a formal realization of this puncture "process", 
where the light-cone becomes degenerate. 

In principle we can now calculate the process where we start out with m 
spatial universes at time and end with n universes at time t, represented 
as 



Gx,gXh, .., U l[, .., C; t) = (Om'i) ■ ■ ■ e-*^*t(Zi) . . . ^^IJIO). (39) 



5.1 Dyson- Schwinger equations 

The disk amphtude is one of a set of functions for which it is possible to derive 
Dyson-Schwingcr equations (DSE). The disk amplitude is characterized by 
the fact that at t = we have a spatial universe of some length, and at some 
point it vanishes in the "vacuum" . Let us consider the more general situation 
where a set of spatial universes of some lengths k exists at time t — and 
where the universes vanish at later times. Define the generating function: 




(40) 



Notice that if the tadpole term had not been present in H Z{J) would triv- 
ially be equal 1. We have 



lim (0|e' 



*t(z,)...^t(g|o) 



5"Z(J) 



(41) 
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Z{J) is the generating functional for universes that disappear in the vacuum. 
We now have 

d 



= hm 

t—*oo 



(0|e-*^ e/d/ J(0*t(0|o) ^ _(o|e-*^ H e/<^'-^«*'«|0) 



(42) 



Commuting the ^(Z)'s in H past the source term effectively replaces these 
operators by IJ{1), after which they can be moved to the left of any \E'^(/) 
and outside (0|. After that the remaining ^'^(/)'s in H can be replaced by 
S/SJ{1) and also moved outside (0|, leaving us with a integro-differential 
operator acting on Z{J): 

= ^ dU{l)o(^l,J,jj^Z{J) (43) 

where 

O(u,^)=«o(i)j^-W (44) 

Z{J) is a generating functional which also includes totally disconnected 
universes which never "interact" with each other. The generating functional 
for connected universes is obtained in the standard way from field theory by 
taking the logarithm of Z[J). Thus we write: 

F(J) = logZ(J), (45) 

and we have 



lim(0|a-'%t(«---**(«|OU=jj(^P^ 



(46) 



J=0 

and we can readily transfer the DSE (43)- (44) into an equation for the con- 
nected functional -F(J): 

,^[di AD [h^H) ^ - sii) - ,j [ dl' jjgj^^ 

, f',,sF{j) sFjj) , r rnr. ^nj) \ ,,7, 

w^m^)' Jo ^^^^^^^5j(7T]or ^ ^ 

Prom eq. (47) one obtains the DSE by differentiating (47) after J(l) a number 
of times and then taking J{1) — 0. 
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5.2 Application of the DSE 

Let us introduce the notation 

5"F(J) 



(48) 



5J(h)---5J{ln) 

as well as the Laplace transform w{xi, . . . , Xn)- Let us differentiate eq. (47) 
after J{1) one and two times, then take J{1) — and Laplace transform the 
obtained equations. We obtain the following equations (where Ho{x)f{x) — 

d^ix' - x)f{x)]y. 

= Ho{x)w{x) — 1 + gsdx(w{x,x) + w{x)w{x)j , (49) 

= {Ho{x) + Ho{y))w{x,y)+gsd^w{x,x,y) + gsdyw{x,y,y) (50) 
+2gs {d^[w{x)w{x, y)]+dy[w{y)w{x, y)]) + 2agsd^dy /' ^(^) 



X — y ) 

The structure of the DSE for an increasing number of arguments is hopefully 
clear (see [16] for details). 

We can solve the DSE iteratively. For this purpose let us introduce the 
expansion of w{xi, . . . , x„) in terms of the coupling constants gs and a: 

oo oo 

w{xi,...,Xn)^ ^ Oi'' ^ g"^ w{xu...,Xn]m,k). (51) 

k=n—l m=k—l 

The amplitude w{xi, . . . , x„) starts with the power {ag,^)^^^ since we have to 
perform n mergings during the time evolution in order to create a connected 
geometry if we begin with n separated spatial loops. Thus one can find the 
lowest order contribution to 'w{xi) from (49), use that to find the lowest order 
contribution to w(xi,X2) from (50), etc. Returning to eq. (49) we can use 
the lowest order expression for w{xi,X2) to find the next order correction to 
'w{xi), etc. 

As mentioned above the amplitude w{xi, . . . ,Xn) starts with the power 
{ags)"'~^ coming from merging the n disconnected spatial universes. The 
rest of the powers of ags will result in a topology change of the resulting, 
connected worldsheet. From an Euclidean point of view it is thus more 
appropriate to reorganize the series as follows 

oo 

w{xi, . . . , a;„) = (as's)""^ ^{agl)''wh{xi, ■■■,Xn) (52) 

^=0 
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Whixi, ...,Xn) = "^giw^xi, . . . ,Xn;n-l+2h+j,n-l + h) (53) 

j=0 

and aim for a topological expansion in ag^, at each order solving for all 
possible baby-universe creations which at some point will vanish into the 
vacuum. Thus Whixi, ■ ■ ■ ,Xn) will be a function of although we do not 
write it explicitly. The DSE allow us to obtain the topological expansion 
iteratively, much the same way we already did as a power expansion in gg. 



6 The matrix model, once again 

Let us consider our N x N Hermitian matrix with the qubic potential (28) 
and define the observable 

W(.„ .... x„). = iV- ((TV (^) . . . (tr ^) ) . (54) 

where the subscript d refers to the fact that the correlator will contain 
disconnected parts. We denote the connected part of the correlator by 
iy(a;i . . . ,Xn)- It is standard matrix model technology to find the matrix 
model DSEs for W{xi . . . ,x„). We refer to [18, 19, 20, 21] for details. One 
obtains precisely the same set of coupled equations as (49)-(50) if we identify: 

and the discussion surrounding the expansion (52) is nothing but the stan- 
dard discussion of the large N expansion of the multi-loop correlators (see 
for instance [21] or the more recent papers [22, 23, 24]). Thus we conclude 
that there is a perturbative agreement between the matrix model (28) and 
the CDT SFT in the sense that perturbatively: 

W{xi,...,Xn)^w{xi,...,Xn). (56) 

In practice the SFT is only defined perturbatively, although in principle we 
have available the string field Hamiltonian. However, we can now use the 
matrix model to extract non-pertubative information. The identification of 
the matrix model and the CDT SFT DSEs were based on (55), but in the SFT 
we are interested in o; = 1, i.e. formally in = 1, in which case the matrix 
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integrals reduce to ordinary integrals. This means that we will consider the 
entire sum over topologies "in one go" : 

(57) 





1 




y dm exp 




Am m^ 1 


. 9s 


V 3 



Z{9A) 

while the observables (54) can be written as 



exp 



Wd{xi, ...,Xn)^ —, TT / dm r r^. (58) 



Z{gs-, A) J (xi - m) • • • {xn - m) 

These integrals should be understood as formal power series in the dimen- 
sionless variable k defined by eq. (18). Any choice of an integration contour 
which makes the integral well defined and reproduces the formal power series 
is a potential nonperturbative definition of these observables. However, dif- 
ferent contours might produce different nonperturbative contributions (i.e. 
which cannot be expanded in powers of t), and there may even be nonper- 
turbative contributions which are not captured by any choice of integration 
contour. As usual in such situations, additional physics input is needed to 
fix these contributions. 

To illustrate the point, let us start by evaluating the partition function 
given in (57). We have to decide on an integration path in the complex plane 
in order to define the integral. One possibility is to take a path along the 
negative axis and then along either the positive or the negative imaginary 
axis. The corresponding integrals are 

Z(^„A) = a/A«:1/^F±(/^-2/3), F±(«:-2/3) = 27re±^'^/«Ai(«:-2/V2-/3)^ (59) 

where Ai denotes the Airy function. Both F± have the same asymptotic 
expansion in k, with positive coefficients. Had we chosen the integration 
path entirely along the imaginary axis we would have obtained (27ri times) 
Ai(K~^/^), but this has an asymptotic expansion in k with coefficients of 
oscillating sign, which is at odds with its interpretation as a probability 
amplitude. In the notation of [25] we have 

F±(z) = TT (Bi{z) ± iAi(^)) , (60) 

from which one deduces immediately that the functions F±{k^'^^^) arc not 
real. However, since Bi(«;~^/^) grows hke es^ for small k while Ai(K~^/^) falls 
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off like e~3^, their imaginary parts are exponentially small in 1/k compared 

to the real part, and therefore do not contribute to the asymptotic expansion 
in K,. An obvious way to define a partition function which is real and shares 
the same asymptotic expansion is by symmetrization, 

^(F+ + F_) = 7rBi. (61) 

The situation parallels the one encountered in the double scaling limit of the 
"old" matrix model [26], and discussed in detail in [27], but is less compli- 
cated. We will return to a discussion of this in Sec. 8. 

Presently, let us collectively denote by F{z) any of the functions F±[z) 
or TTBi{z), leading to the tentative identification 

Z{g,,X) ^VXk'/'f[k~'/'), F"{z) = zF{z), (62) 

where we have included the difi^erential equation satisfied by the Airy func- 
tions for later reference. Remarkably, this partition function was also found 
in [28] , where a double scahng limit of so-called branched polymers were stud- 
ied. It reflects that a significant part of the dynamics associated with the 
branching, as reflected in flgs. 3 and 4, are indeed captured by a branched 
polymer model. It should be no surprise that this is possible. Branched 
polymers play an important role in non-critical string theory [29, 30] and 
even in higher dimensional Euclidean quantum gravity [31]. However, the 
2d CDT has a much richer set of observables than those encountered in the 
theory of branched polymers, for instance the observables Wd{xi, . . . ,x„), 
the calculation of which we now turn to. 

Let us introduce the dimensionless variables 

x = x\/X, m = gy^(3, Wd{xi, ...,Xn) = X'"^'^Wd{xi, . . . , x„). (63) 

Assuming Xk > 0, we can write 
1 /"°° 

i_^i/sf3 - I d« exp [- {x - «:V3^) . (54) 

We can use this identity to re-express the pole terms in eq. (58) to obtain 
the integral representation 

Wd{x,, . . . , x„) = / rr da, e-(*i"i+-+*"'^") ^ ^ (65) 

•^0 t\ f(^k-I) 
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for the amplitude with dimensionless arguments. By an inverse Laplace 
transformation we thus obtain: 

For the special case n = 1 we find 

(ar) 

for the disc amplitude, together with the remarkable relation 

W,{k,...,ln)=W{k + --- + Q. (68) 
By Laplace transformation this formula implies the relation 

From Wd{xi, . . . ,Xn) we can construct the connected multiloop functions 
W{xi, . . . ,Xn) using standard field theory. Let us remark that the asymp- 
totic expansion in k of W{xi, . . . , Xn) of course agrees with that obtained by 
recursively solving the CDT Dyson- Schwinger equations. 



7 Stochastic quantization 

It is a most remarkable fact that the above mentioned results can all be 
understood as a result of stochastic quantization of space. In this picture 
time becomes the stochastic time related with the branching of space into 
baby universes and the original CDT model described in Sec. 2 becomes the 
classical limit where no stochastic processes are present [32]. 

Recall the Langevin stochastic differential equation for a single variable 
X (see, for example, [33, 34]). 

x(''\t)^-f(x^''\t)^+Vfliy{t), (70) 

where the dot denotes differentiation with respect to stochastic time t, i/{t) is 
a Gaussian white- noise term of unit width and f{x) a dissipative drift force: 
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The noise term creates a probability distribution of x{t), reflecting the as- 
sumed stochastic nature of the noise term, with an associated probabihty 
distribution 

P{x,xo;t) = {6{x-x'^''\t;xo)))^, (72) 

where the expectation value refers to an average over the Gaussian noise. 
P{x, xq; t) satisfies the Fokker-Planck equation 

^ ^ ypiY^-^t) ^ ^^^^^^^^ ^ A _ ^^3^ 



dt dx \2 dx 

This is an imaginary-time Schrodinger equation, with playing a role 
similar to h. It enables us to write P as a propagator for a Hamiltonian 
operator H, 

P{x,xo;t) = (x|e-*^|a;o), H = ^Qp' + tp f{x), (74) 

with initial condition x{t = 0) = xq, and p — —idx- It follows that by 
defining 

d 

G{xo, x; t) = -— P{x, Xq; t) (75) 

OXo 

the function G{xq, x; t) satisfies the differential equation 

dG{xo,x;t) d (l dG{xo,x;t) . 

dt = 1^2 • • ^''^ 

Omitting the noise term corresponds to taking the limit Q — > 0. One can 
then drop the functional average over the noise in (72) to obtain 

d 

Pci{x,xo;t) ^ S{x - x{t,xo)), Gci{xo,x;t) ^ ——5{x - x{t,xo)). (77) 

OXo 

It is readily verified that these functions satisfy eqs. (73) and (76) with f2 = 0. 
Thus we have for S{x) — —Xx + x^/3: 



dGci{xo,x;t) d 



(^{X-xl)Ga{xo,x;t)y (78) 



dt dxo 

Comparing now eqs. (7) and (78), we see that wc can formally re- interpret 
Gf^ {xq, x;t) - an amplitude obtained by nonperturbatively quantizing Lorentzian 
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pure gravity in two dimensions - as the ''classical probability" Gci{xo,x;t) 
corresponding to the action S{x) = ~Xx + x^ /3 of a zero-dimensional system 
in the context of stochastic quantization, only is the boundary condition dif- 
ferent, since in the case of CDT x is not an ordinary real variable, but the 
cosmological constant. The correct boundary conditions are thus the ones 
stated in eqs. (8)-(9) 

Stochastic quantization of the system amounts to replacing 

Gf\xo,x;t)^Gixo,x;t), (79) 

where G{xo,x;t) satisfies the differential equation corresponding to eq. (76), 
namely. 

We have introduced the parameter gg :— ^2/2, which will allow us to repro- 
duce the matrix model and SFT results reported above. 

A neat geometric interpretation of how stochastic quantization can cap- 
ture topologically nontrivial amplitudes has been given in [35]. Applied to 
the present case, we can view the propagation in stochastic time t for a given 
noise term i/{t) as classical in the sense that solving the Langevin equation 
(70) for x'^^'^t) iteratively gives precisely the tree diagrams with one exter- 
nal leg corresponding to the action S{x) (and including the derivative term 
x^'^\t)), with the noise term acting as a source term. Performing the func- 
tional integration over the Gaussian noise term corresponds to integrating 
out the sources and creating loops, or, if we have several independent trees, 
to merging these trees and creating diagrams with several external legs. If 
the dynamics of the quantum states of the spatial universe takes place via the 
strictly causal CDT-propagator Gq = e~^^° , a single spatial universe of length 
/ cannot split into two spatial universes. Similarly, no two spatial universes 
are allowed to merge as a function of stochastic time. However, introduc- 
ing the noise term and subsequently performing a functional integration over 
it makes these processes possible. This explains how the stochastic quanti- 
zation can automatically generate the amplitudes which are introduced by 
hand in a string field theory, be it of Euclidean character as described in [35] , 
or within the framework of CDT. 

What is new in the CDT string field theory considered here is that we 
can use the corresponding stochastic field theory to solve the model, since 
we arrive at closed equations valid to all orders in the genus expansion. Let 
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us translate equations (80) to Z-space 



dG{lo,l;t) 
dt 



H{lo)G{lo,l;t), 



(81) 



where the extended Hamiltonian 



H{1) 



(82) 



now has an extra potential term coming from the inclusion of branching 
points compared to the Hamiltonian Ho(l) defined in (14). It is truely re- 
markable that all branching and joining is contained in this simple extra 
term. Formally H{1) is a well-defined Hermitian operator with respect to the 
measure (32) (we will discuss some subtleties in the next section). 

We can now write down the generalization of Wheeler- de Witt equation 
(32) for the disk amplitude 



Contrary to W^^^l) appearing in (32), W{1) contains all branchings and all 
topology changes, and the solution is precisely (67)! This justifies the choice 
Qs = fl/2 mentioned above. Recall that E = does not belong to the 
spectrum of Ho{l) since Wq{1) is not integrable at zero with respect to the 
measure (32). Exactly the same is true for the extended Hamiltonian H{1) 
and the corresponding Hartle-Hawking amplitude W{1) 
We have also as a generahzation of (13) that 



describes the nonperturbative propagation of a spatial loop of length Iq to a 
spatial loop of length / in proper (or stochastic) time t, now including the 
summation over all genera. 

8 The extended Hamiltonian 

In order to analyze the spectrum of H{1), it is convenient to put the differ- 
ential operator into standard form. After a change of variables 



H{l)W{l) = 0. 



(83) 



G{lo,l;t) = {l\e 




(84) 




(85) 
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the eigenvalue equation becomes 

H{z)<t>{z) = E<p{z), H{z) = + \\z' + A _ |^4_ ^gg) 

This shows that the potential is unbounded from below, but such that the 
eigenvalue spectrum is still discrete: whenever the potential is unbounded 
below with fall-off faster than —2;^, the spectrum is discrete, reflecting the 
fact that the classical escape time to infinity is finite (see [36] for a detailed 
discussion relevant to the present situation). For small gs, there is a large 
barrier of height A^/ {2gs) separating the unboTinded region for / > \/ Qg from 
the region < / < \/{2gs) where the potential grows. This situation is 
perfectly suited to applying a standard WKB analysis. For energies less 
than \^/{2gs)i the eigenfunctions (33) of Hq{1) will be good approximations 
to those of H{1). However, when I > \/gs the exponential fall-off of i^n\l) 
will be replaced by an oscillatory behaviour, with the wave function falling 
off only like The corresponding ipn{l) is still square-integrable since we 

have to use the measure (32). For energies larger than A^/ {2gs)i the solutions 
will be entirely oscillatory, but still square-integrable. 

Thus a somewhat drastic change has occurred in the quantum behaviour 
of the one-dimensional universe as a consequence of allowing topology changes. 
In the original, strictly causal quantum gravity model an eigenstate il)n\l) 
of the spatial universe had an average size of order 1 / \/~\. However, allowing 
for branching and topology change, the average size of the universe is now 
infinite! 

As discussed in [36], Hamiltonians with unbounded potentials like (86) 
have a one-parameter family of selfadjoint extensions and we still have to 
choose one of those such that the spectrum of H{1) can be determined un- 
ambiguously. One way of doing this is to appeal again to stochastic quan- 
tization, following the strategy used by Greensite and Halpern [37], which 
was applied to the double-scaling hmit of matrix models in [38, 39, 36]. The 
Hamiltonian (74) corresponding to the Fokker-Planck equation (80), namely, 

i/(.)^(x) = -,.^ + -(^^^(.)j, Six)=[^--Xxym 

is not Hermitian if we view x as an ordinary real variable and wave functions 
'4){x) as endowed with the standard scalar product on the real line. However, 
by a similarity transformation one can transform H{x) to a new operator 
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which is Hermitian on L'^{R, dx). We have 

which after substitution of the exphcit form of the action becomes 

H{^) = -Qs^ + V{x), V{x) = ^(A - x^f + X. (90) 
The fact that one can write 

«.-v^iL + ^l^ (91) 

implies that the spectrum of H{x) is positive, discrete and unambiguous. We 
conclude that the formalism of stochastic quantization has provided us with 
a nonperturbative definition of the CDT string field theory. 
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